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ABSTRACT 


An analytic model is developed to allow derivation of the 
probability that a molecule in the atmosphere of Saturn's 
rings collides with at least one ring particle when tra- 
versing the ring plane. The resulting expression involves 
details of the molecule's trajectory, including the velocity 
of the molecule relative to ring particles; thus, the theory 
is compatible with recently developed ballistic-transport 
computer models used in the analysis of the rings. The 
Souelrsion theory is applied to the case of a low energy, 
isotropic molecular production process to make inferences on 
the spatial extent of the ring atmosphere associated with 
such a source. The high frequency of collisions expected 
for the atmosphere in the vicinity of the A and B rings 


suggests a toroidal atmosphere. 
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LIST OF “Ss 2SOrs 


thickness of a sublayer of a "slab" whose overall 
thickness is h (h < b) 


thickness of the translating slab (containing par- 
ticles) through which a molecule passes (corres- 
ponds to ring thickness in application oreey 
collision theory) 


number density of particles ina slab (or in the 
ring for applications of the collision theomg 


particle radius (sometimes referred to in this 
paper as the "size" of a particle) 


ejection velocity of a molecule, relative to its 
source particle 


velocity of a molecule 


velocity with which the slab translates (equivalent 
to a representative) velocity of the particles 


"“COntained” in the slab 
distance that the slab translates in time t 


scale length for the extent of the slab in direc- 
tions normal to its thickness 


equilibrium number of molecules (of a given type) 
in ring atmosphere deriving from a specific produc} 
tion mechanism 


number of slab (ring) particles of size £5) eom— 
tained in the slab (meaningful only for discretely 
distributed r) 


probability that a molecule collides with at least 
one ring particle during a given passage through 
the slab (ring) 


probability that a molecule passes through the 
slab (ring) without colliding with a parervelomas 
size cr; (only strictly meaningful for discerecoum 


probability that j particles of size r. are not in 
the collisien region 


@e(r,A) 


extinction cross section, sometimes referred to as 
the extinction coefficient (relates to the effec- 

tive size of a particle with regard to removal of 

electromagnetic energy from a beam via scattering 

and absorption) 


orbital radius (see Section III) 
interaction radius (see Appendix B) 


sticking coefficent (the probability that a mole- 
cule will adhere to a surface when it collides 
with that surface) 


a "mean" orbital period typical of the trajectories 
of all molecules of a certain type produced by a 
given process 


number of completed orbital periods 


production rate (of molecules) associated with a 
specific mechanism, usually including only those 
molecules which assume closed orbits 


direction cosines which specify a molecule'’s path 
as it passes through the slab, when working in 

a frame of reference in which the slab translates 
at ws (SS 2 KC Ue ae) 


direction cosines which specify a molecule's path 
through the slab, when working in a frame of refer- 
ence in which the slab is fixed (see Figure C2) 


aiicgiiomused = cordescribe the deviation of the actual 
direction of v_ from the assumption that v_ is 
COnseant p 


wavelength (of electromagnetic radiation) 


angle that an electromagnetic beam makes with the 
Sumbace: Memmal Of a slab (ring) 


BcCie 1 Heemoee=y plane between the projection of 
the molecular path and the +x direction (v 
direction), defined in a frame of reference in 
which v, # 0, referred to in this paper as the 
"orbital aximuth" 


same as 8, except defined in the frame of reference 
with Me = 0 
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I. SN ERODGET ION 


It has been well established that the particles of which 
Saturn's rings are composed are primarily water ice (Cuzzi 
et al., 1984). Accepting this, one would expect that various 
erosive processes would produce a ring atmosphere comprised 
hig HO molecules and their constituents, e.g., H atoms. The 
density of such an atmosphere would be determined by balancing 
Mrempreduction rates {of atoms and molecules) associated with 
the erosion mechanisms with losses resulting from capture by 
Sewn, eScCape from the Saturnian system, and recapture by the 
MEE, particles. 

Farly theoretical estimates of a ring atmosphere (Denne- 
feld, 1974; Blamont, 1974) considered sources such as subli- 
mation, meteroid impact, and bombardment by solar and 
interstellar wind.+ These analyses were concerned only with 
the atmosphere associated with the A and B rings (see Figure l), 
by far the most significant potential sources of material for 
a ring atmosphere, owing to their large masses (Appendix A 
gives a summary of ring nomenclature and dimensions). Denne- 
Feld calculated the density of the H 


2 
em, the molecules confined to a toroidal region closely 


O atmosphere to be ~25 


surrounding the A and B rings. He estimated the H atom 


“Dennefeld presented a conceptual model for the study of 
ring atmospheres, which is of fundamental importance. This 
model is discussed later in this section. 
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Figure l. 


Nomenclature of the Sadturnian Ring s7oeem 
(drawn to scale) 


component of the atmosphere to have a density of ~l cm > 


distributed tChreugh a. spnerer ons x 10> km, well beyond the 
outside edge of the A ring. 


Schemes for detection of the atmosphere all rely on 


measurement of H Lyman-a emission from the vicinity of the 
rings, the most sensitive means available. Thus, analysis is 
limited to the neutral H atmosphere. Using this approach, 


several researchers put Dennefeld's estimates to the test. 


Observations from instruments onboard rockets (Weiser et al., 
1977), satellites 


(Barker, et al., 1980; Clarke eo eee 


LO 


Puoneea eth (Judge et als, 1980), and Voyager I (Broadfoot et 
al., 1981) detected H atoms near the A and B rings. Results 
were reasonably consistent, Suggesting a H density 2-3 orders 
of magnitude greater than Dennefeld predicted. The shape 

of the atmosphere (spherical versus toroidal) cannot be con- 
clusively deduced from these Lyman-a measurements, but pre- 
liminary analysis of Voyager data indicates that it does not 
extend beyond the outer edge of the A ring (Cuzzi et al., 
moe) . 

Efforts to reconcile theoretical estimates of the H 
atmosphere with observations have primarily been directed 
toward proposing new sources or re-evaluating those previously 
considered. Proton flux from Saturn's ionosphere (Ip, 1978), 
photodissociation of H,Q on the ring particles (Carlson, 1980), 
Magnetospheric ion sputtering (Cheng et al., 1978,1980), and 
Meteroid impact (Morfill et al., 1983; Ip, 1983) are among 
the most seriously studied erosive mechanisms. Of these 
meteroid impact appears likely to be the dominant source. 

While understanding of the ring atmosphere has increased, 
Dennefeld's conceptual model has remained the foundation for 
its theoretical analysis. It is worthwhile to review this 
model. Suppose that a particular erosive process produces 


atmospheric molecules (e.g., H50 or H) with these molecules 


“some have also considered H atom sources external to the 
rings such as atoms escaping from Saturn's atmosphere 
(Shemansky and Smith, 1982). Preliminary calculations indi- 
cate that this process may be quite significant. 


La 


having some initial velocity distribution. @he trajeceoms, 
of each molecule is primarily governed by its gravitational 
interaction with Saturn, the dominant central Sedyenieee 
Saturnian system. If the ring atmosphere is quite tenuous, 
collisions between molecules as they move in their Keplerian 
Orbits is unlikely (this assumption poses no problem for 
Saturn's rings, since for a density of even 1000 em the 
mean free path is 107 km). A molecule's initial velocity 
and position within the ring determines its orbit for which 
there are three possible cases* 


(1) the trajectory takes the molecule into Saturn's 
atmosphere, thus removing it from the ring atmosphere. 


(2).. the orbit is open (hyperbolic) with the molecule 
escaping the entire Saturnian system. 


(3) the orbit is elliptic with the molecule passing 
through the ring once or twice during each orbital 
period (depending on the eccentricity of the orbit). 

For those molecules with elliptic orbits, some will be lost 
from the atmosphere to the rings. The probability that this 
occurs during each passage of a molecule through the rings 
1s the product of the probabilities of its colliding wie 
ring particles and of its sticking to a particle during wom. 
process. The spatial extent of the atmosphere tends to be 


toroidal if the initial speed of the molecules relative to 


the orbital speed of the ring particles is small. On the 


3The first two cases apply if the erosive mechanism 
produces molecules with high speeds relative to the ring 
particles; the third is dominant for relatively lower energy 
processes. 
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other hand if the initial speed is much greater than that of 
the ring particles, then the atmosphere will tend to be more 
spherical. A mathematical expression of this process of the 


form introduced by Dennefeld is 


where: 


N number of molecules (of a given "kind") in 
the equilibrium atmosphere (from a specific 


production mechanism). 


W = production rate (sec a of molecules asso- 
Ciated with a specific mechanism which 
7 assume elliptic orbits (call this the 
"source term"). 


T = mean orbital period (sec); this term incor- 
Bomacteseimiieta| velocity information, so it 
might more generally be referred to as the 
Dae ist rc. tema 


Seco ecomng acocthicientjethe probability that 
a molecule will adhere to a ring particle 
when it collides with it. 


P = probability that the molecule collides with 
a eee particle as it passes through the 
iE sWgvel 


This model is admittedly crude, but it does contain the 
essence of the physical system. To see how the theory has 


evolved, each term is discussed below. 


4 tacit assumption in this formulation is that the 
average trajectory intersects the ring only once during 
an orbital period. If it passes through the ring twice 
@@eing a period, P, hasa slightly different interpretation, 
but the other terms are unaffected. 


IES: 


The source term, W, has received much attention as dis- 
cussed in some detail earlier. 

The ballistic term, os receiving increased emphasis. 
Ballistic transport models (Durisen, 1984; Ip, 1983) aneegme 
porate inhomogeneities found in the macrostructure of the 
rings with improved estimates of initial velocity distribu- 
tions to allow better specification of molecular trajectories. 

The sticking coefficient, S, is probably the least under- 
stood factor in the theory. The physical chemistry under- 
lying the phenomenon is not well developed, while laboratory 
experiments: do not appear to accurately replicate conditions 
Found -1n the ring environment.> While it may be possible that 
the physical chemistry problem will be solved, a more basic 
dilemma will likely remain for quite some time: S cannot be 
determined accurately unless detailed knowledge of individual 
ring particles is developed. Specifically, impurity levels 
and the phase of the water ice in the particles are poorly 
known (Weidenschilling, 1984), and the surface texture of 
the particles is still being debated (Kerr, 1985). 

Another feature of the model, which has received little 


attention is the collision prebaliliez i In the literature 


>For example, Carlson (1980) used §S = 0.22 for collisions 
involving H atoms, based on experiments (Brackmann and Fite, 
1961) in which a H beam impinged upon a liquid-nitrogen— 
cooled copper plate. The beam was so dense that recombination 
of atoms to produce Hz was a significant source of loss OF 
the atoms from the beam. This experiment does not appear 
to closely approximate the physical system. 
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Ey is consistently expressed as a function of the optical 
thickness, tT, the definition of which is discussed in Section 
III.B. Dennefeld set Ee equal to tT, while recent works 


fearison, 1980; Ip, 1983) use 


This seems to be a simplification of the collision process, 
Mee one would expect es to depend on such trajectory-related 
factors as velocity of the molecule relative to the ring 
particles at the location of passage, and the path length 
through the ring (Durisen (1984) states that ee will depend 
On optical thickness, and the "slant path of an ejectum 
through the rings," but does not elaborate). Also, the above 
expressions for es do not account for multiple collisions 
during one passage through the rings by a molecule. ° The 
existence of multiple collisions, if significantly probable, 
has implications both for density of the atmosphere (there 
are more "sticking" opportunities), and for its spatial 
m@estribution. 

The aim of this paper is to elevate the theory of colli- 


Sions between ring particles and atmospheric molecules to a 


Me7el of sophistication on par with that of ballistic 


eNote that nowhere in Dennefeld's model are multiple 


collisions considered. The impact of this omission is 
discussed later. 


i> 


transport models now utilized in the study of ring atmos- 
pheres. A generalized model for determining Ee is developed 
in Section II, while in Section III its applveas te ae 

the actual ring environment is discussed. Section IV is a 
mathematical analysis of the expression derived for aoe 
Finally, in Section V the theory is applied to the speeuem 
case of low energy, isotropic production processes (e.g., 
sublimation) to illustrate inferences that can be made on 
the spatial extent of the ring atmosphere, based on the fre- 


quency of collisions. 


as 


Pie pee aeORy 


In this section, the general problem of a molecule passing 
through a translating target consisting of randomly dis- 
tributed macroscopic particles is addressed. The proba- 
bility that it collides with at least one particle is 
earculated. | 

Consider a region in space of large extent, L, in two 
dimensions, and of a smaller thickness, h. This slab is 
populated by noninteracting spherical particles which are 
randomly distributed (spatially) and fixed with respect 
to the slab. Assume that the size of the spheres is limited 
to radii in the interval Loan! ic << h] upon which the 
radii may be distributed either discretely or continuously. 
The number density of particles of size r is denoted by 
me). Finally, the slab as a whole translates in a direction 
normal to its thickness with constant velocity, ve 

Now suppose that a molecule traveling along a straight 
path with speed Vv, passes through the slab. Assume that the 
molecule is much smaller than the smallest particle in the 
Slab, and that it only interacts with the particles through 


mechanical collision. The relationship between the molecule 


Tor discrete distributions n(r) has units m > For 
Continuous distributions it has units m74, so that n(r)dr 
iS interpreted as the number of particles per m’ in a radius 
increment [r,rt+dr]. 


Ly 


and the slab can be conveniently described using direcourem 
cosines as shown in Figure 2(a). Notice that the x-axis is 


chosen to be parallel to v 





L 


(a) v, # 0 (b ) vist 


Figure 2. Alternate Frames of Reference 


Calculation of the probability of the molecule colliding 
with at least one particle as it passes through the slab, 
henceforth denoted as Pas can be accomplished in the tram 
of reference shown in Figure 2(a). This procedure, however, 
is mathematically tedious (see Appendix B). The problem is 
Simplified by converting to a frame of reference in which 
We equals zero (Figure 2(b)). This approach is Getailed 
below. 

Suppose that the conversion to the frame Cf rererenee 


with ws equal to zero has already been accomplished (see 


i 


eeeecnaix C for details). For simplicity, let the path of 

the molecule be normal to the slab, i.e., y' equals zero, and 
n(r) be discrete (extension of the results for this case will 
be generalized to cases with continuous r and arbitrary 
molecular path). The molecule will collide with any parti- 
cle of radius re which has its center anywhere in the cylin- 
drical volume (collision region) whose axis coincides with 


the molecule's path as shown in Figure as The probability 


ee molecular path (normal to 
1 Ene Ts Faby 


collision region 





Figure 3. Collision Region for Particles of Radius r 


meaat a COllision does not occur with particles of size cr. 


memequivalent to the probability that no particles of that 


SNote that it was tacitly assumed that the molecule is 
essentially a point mass. This enables one to define the 
collision region as was done above. An equivalent way to 
Meee ptualize the Collision cylinder is to let the mole- 
cule have radius r, and the particles be point masses. The 
moebecule then collides with any particles it "sweeps out" 
along its path. Of course the resulting collision region 
again 1s a cylinder of radius r. 


IE, 


Size have centers in the collision region. Denote this 
probability as Pio (ky): Since the particles were assumed to 


NOt win teraceaaele Pg (¥4? for all i are mutually independent. 


This allows one to write £0r diserecreun. 


P = l-MIP 
a 


é Oo (hy). over all rs (lik 


gl 


The problem is to find an expression for Pig (ty) - 

Under the assumptions of the model, the uniform proba- 
bility distribution governs the location of each partieue 
1.e., the probability that a particle will be found in some 
portion of the slab is proportional to the volume of that 


portion. For example, consider one particle of size r.. 


The probability that it is not in the collision region 


where the numerator is volume of the slab minus the collision 
region, and the denominator is the volume of the slab. Since 
the location of each particle in the slab is independéeneaak 
the location of other particles, the probability thacewaag 

two specific particles of size r, are not in the collisvem 


region 1s 


tg 
li 
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Ppeeending this line Of reasOnming, it is clear that 


é 2 
Lh -aryh N(r,) 


a 


Po. bree = 
= L“h 


where: 


is the total number of particles with radius rs in the slab. 


fmeeeir.) can be rewritten as 
mie CL 


mein (x, )h N(x.) 


Lae aa = (ae = See katiee2)) 


ne N(r,) is large, then 
_ 2 
Peer) = exp{-mr;n(r,)h} : aes) 


SMestituting into equation (II.1) with N(r,) large for all 


rc; yields 
P el exp{-mrsn(r a) 
c ; 1 a 
one 
P = 1 - expt-th Yren(r ei (II.4) 
@ 7 1 1 : ; 


ZA 


recalling that r is discretely distributed, and the molecuwam 
path is normal. 

Extending this result to the case when the molecule's 
path makes an angle y' with a normal from the slab is easily 
accomplished. Notice that the only difference between this 
case and that of the normal path is that the length of the 


molecular path (hence the length of the collision cylinder) 


increases from h to h/|cos y'|. Therefore, 
Pe = il ©XP\TSos y'] Ma ha (I Ts 


with-r discrete and the molecular path. not necessarily normal. 
Extension to the case when r is continuous is affected 

by replacing the summation in equation (II.5) with integra- 

tion (see the first footnote on page 9 to see the new inter- 


pretation required for n(r)). This yields 


aX 5 
ented ia (II.6) 


=-—7h “m 
Po = 1 = ‘expi — f 
Cc feos y'| z 
min 
where r is continuous and the molecular path is not neces- 
sarily normal to the slab. 
Equations (I1I.5) and (11.6) represent solutions t€0 he 

problem when working ina frame of reference for which V5 


equals zero. Often problems are not initially formulated 


in this particular frame. Appendix C shows how to convert 


ae 


these results back to the original frame of reference (shown 


moerigure 1). Substituting from equation (C.8) yields 
P = 1] -exptf allie Pees. (REN Yr-n(r.)} (lee) 
Cc [cos >| 7 ve : Hl i ; 


for r discrete, and 


= 5 we D max 5 
xe = | —&*Pt ices 1] ag ama i J ig igulce isbed 
iq vif 
Gio} 
meer COMtinuous. Equations (II.7) and (II.8) are more con- 


meatent than equations (II.5) and (11.6) for some applications. 
It is interesting to note that while the length of the 
molecular path in the frame of reference with vo not equal 
to zero is h/|cos y|, the “effective path length," pertinent 


@eethe Collision problem, is 


= = YP YP) 2 
Loge eoEsT IL Pig | ae a + ( ae ’ Ci. 9') 
which accounts for the relative motion of the molecule versus 
mgieetring particles. 
Equations (II.5)-(II.8) represent the fruition of this 
@erivation of Poe Table 1 gives a review of the tacit and 


explicit assumptions of the model. 
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TABE Eu 


ASSUMPTIONS OF THE THEORETICAL e2as 


Assumptions 
spherical particles 


random spatial distribu- 
tions of the particles 


Jee Saat i 


Bee > > ae) a 
5 molecule 


molecule enters and leaves 
slab emo ER its large 
(area of L*) faces 


Slab translates with 
velocity, Set Horna. £oO 
its thickness, h 
v and Vv constant 
p m 
particles fixed relative 
to the slab 


number of particles in the 
slab is large 


Comments 
mathematically convenient 


justifies using the uniform 
probability density fune@men 


facilitates application oes 
the model to the rings 


no Significance beyond 
"quantifying" particle Siig 


effectively excludes colli- 
sions outside the nominal 
confines of the slab 


allows the molecule to be 
treated as a point particle 


facilitates application "es 


the theory to the rings 


facilitates application of 
the theory to the rings, 
mathematically convenient 
mathematically convenient 
justifies exponential 


approximation in probabiGiie 
calculat rem 


Is the theoretical model utilized here relevant to condi= 
tions in Saturn's rings? This question is discussed ie. 


Nex=  Sececelon. 
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Tite) SePEEGATION TO SATURN'S RINGS 


Pees ePLICABILITY OF THE THEORETICAL MODEL 


ili 


General Comments 


Before proceeding with analysis of collisions between 


ring particles and atmospheric molecules, the validity of the 


generalized "slab" model of Section II must be established 


‘ 


vis-a- 


felons” 


i) 


(2}- 


Vis the actual ring environment. The essential assump- 
of the model are: 
the particles in the slab are spherical. 


the radius (size) of any particle in the slab is 
in an interval [Ypin-Tmaxi where rmi, is much 
Geeteeo- ehanecne Size or ene molecule and rnay 
is much smaller than the slab thickness, Nh. 


the spatial distribution of particles within the 

Slab is random, i.e., the location of each particle 

is governed by the uniform probability density function 
with all particles mutually independent (in a sta- 
tistical sense). 


the number of particles within the slab is large 
enough for the exponential approximation (equation 
tiie) tO be valid. 

the velocity of the slab, vo) is constant. 


the velocity of the molecule, Vin? Tomeoms came - 


the slab particles are fixed relative to the motion 
of the slab as a whole. 


” Some of the assumptions listed in Table 1 were adopted 
to facilitate application of the model to the rings. These 
assumptions have no physical significance, so they are not 
discussed here. 


yee) 


If a “physical slab" can be defined in the rangsewaren 
reasonably satisfies the above conditions, then the appli- 
cability of the model is verified. Note, however, that the 
physical slab is an abstract concept, since it is related wo 
the passage of a specific molecule, on a specific trajecromm 
To facilitate the discussion of the above seven 
assumptions, they can be divided into three categories, those 


related to: 


(1) particle properties (assumptions (1) and (2)). 

(2) particle distribution properties (assumptions (3) 
and (4)). 

(3) velocity-related properties (assumptions (5)-(7)). 


Beiioumeacn of these categories is discussed with respect to 
Saturn's A and B rings. The other rings will not be con- 
Sidered, since they are not expected to contribute much to 
the ring atmosphere (owing to their low masses). Observational 
data referenced in Section I supports this exclusion. 
2. Particle Properties 

Individual ring particles have not vet been observed, 
thus absolute verification of the assumption of spherical 
shape cannot be accomplished. However, such a shape is 
adopted ubiquitously in the development of ring models. 
Observational data does not suggest that this traditional 
approach is invalid (Weidenschilling et al., 1984), in fact 
the Voyager radio occultation experiment (Tyler et al., 1983) 


eliminates the possibility of very elongated particles. 
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Therefore, at present the most reasonable shape to attribute 
Bemrcing particles is spherical. 

Most ring particles have sizes in a continuous radius 
[mamerval (Cuzzi et al., 1984; Esposito et al., 1984; 


Beeecnschilling et al.,; 1984) 


Clearly, ee is much greater than the size of a H50 Onaiel 
molecule. The thickness, h, of the A and B rings is nominally 
not larger than 200 m, but locally may be as small as 10 m 
(Cuzzi et al., 1984). The requirement that eae be much less 
meat h should not create a problem for application of the 
@ollision theory. 
3. Spatial Distribution of the Ring Particles 

The simplest way to describe the rings is as a 
collection of mutually colliding particles whose trajectories 
are dominated by the gravitational influence of Saturn. Under 
such conditions the notion of a completely random particle 
distribution throughout the rings is reasonable. The physi- 
cal slab could then be defined as any translating region 
containing the molecular path, provided that the velocity- 
related assumptions are valid for that specific region. 
Meese tne largest such region consistent with the other 
assumptions to be the physical slab. As will be shown in 


the discussion of velocity-related requirements, this region 


AG 


will have a considerable volume. This is desirable, for even 
if n(r) is quite small (as limited measurements have indi- 
cated (Marouf et al., 1983)), the requirement that the ni@mieen 
of particles in the slab be large is likely to be satisfied. 

In actuality the rings exhibit considerable fine 
structure (1.e., variations in optical thickness) on Traduza 
scales down to ~l km (Esposito et al., 1984). Apparenae 
the primary cause of this fine structure is linked to the 
gravitational effects of other bodies in the Saturnian system, 
especially the moons. Gravitational resonances are respon- 
sible for the creation of zones where stable particle orbits 
are precluded (e.g., the Cassini Division), and may Genememae 
density waves within the rings. The net effect of fine struc- 
ture on the collision theory is to limit the scale over which 
random spatial particle distribution is correct. This cCangee 
best understood by considering the following examples. 

In the first example, the slab is divided into two 


regions as shown in Figure 4. Suppose that the particles in 


molecular 










path 
Note: For simpliciegs 
consider the franemem 
reference with v.. = 0, 

h and require the folecular 


path to be paralleliige 
a slab edge. 





K— 5 Lb —*# —_. 5. —4 


Figure 4. Model for AnalyZing the Effect of Nonrangem 
Particle Distmimemri1on 
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each region are constrained to remain within that region; 

this represents the departure from truly random spatial dis- 
frabution of the particles. Let n, (rx) equal n,(r) with 

r discrete. What is the probability that the molecule 

suffers at least one collision when passing through the entire 
Slab? Attack this problem by reformulating it as passage 


through two smaller slabs shown in Figure 5 as shaded regions. 





Figure 5. The Two Slabs used in Reformulating the Problem 


meem section Ii, the probability of no collision with parti- 


cles of size rc. fOr the Slab of thickness b is 


where: 


For ail re this becomes 


5 ees aCe ) Tess TT, N(x. 


no (oy : N(r, 


which for large N(x.) becomes 


_ 2 b 
Pno(b) 7 7 exptoani niga’ {cou gee 


-Tb 2 
exP (eos all i ren(x,)} . (l7Ti aa 


Similarly, for the slabvorf “echicknmess ) (ia) 


tT (h-b) 


aGerr L ig = oheee ah (ITT aay 


Oia ae oe 


fOr n(x.) L* (h-b) large. The probability of collision fOr 


entire passage through the original slab is 


E (Ill See 


C ~ Pao (b) *no(h-b) 


Substitution yields 


dese) 


Pp =il1- comets ) r“n(xr.)} } r°n(r.) } 
el ae lcos y'] ; 1 ae expt cos y a 1a : 
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Ox. 


= =a 2 
ae = l1- e=P TT cosmy| } r.n(x,)} : le he) 


This result is the same as that obtained when a random spatial 
Pict ibution over the entire slab is assumed. Careful con- 
sideration reveals an important qualification to this conclu- 
sion. Notice that the assumption of exponential form is less 
tenable in the "nonrandom" case, Since there it 1s reguired 
that both n(r, \L“b and n(x, ee oem nile an che 
completely random case en must be large. Extending this 
line of reasoning to cases where the original slab is further 
subdivided into regions to which the particles are confined, 
it is clear that at some point the exponential approximation 
feel certainly fail. 

In the second sane e, the slab is again divided into 
two regions (Figure 4), but now n, (4) is not equal to n5(r). 


Following the approach taken above, the collision probability 


is 


(h-b 
a ~ewtrese yy ) i eewtgs yt) Ps (rj), 


Ong 


_ s =a) Z 
Po = ] BN eee a te } ren, (r,) 


Care 


yes ; (alice 5) 


+ (h-b) } ren 
5 J 


fal 


peed tie) > = 
——— 


where both ees ue and ny (r5)L° (h-b) must be large enough 
to justify the exponential approximation. The collisitea 
theory is here valid in principle, but its application as 
MOGerair hel te 

From this point on, the A and B rings will be assumed 
to satisfy the spatial distribution conditions. Howevergume 
must be kept in mind that fine structure effects may inter- 
Jece inaccunraciesmanto calemiaerens. 

4. Velocity-related Assumptions 

Reviewing, the three velocity-related assumptions 
are: v5 is constant; a is constant; and the particles in the 
slab are fixed relative to vo: Each in turn is addressed 
below, but the validity of the concept of slab velocity when 
applied to the ring system must first be demonstrated. As 
discussed in Section III.A.3, choosing the physical slalwas 
be as large as possible is advantageous. By means of a simple 
example, it can be shown that a single velocity well approxi- 
mates the orbital motion of all of the particles in even a 
quite large slab (this velocity corresponds to v,)- Suppose 
that a slab with L equal to 10,000 km is centered in Enoms 
ring as shown in Figure 6. If us is set equal to the parti- 
cle velocity at C, Vo how much error is inherent to the 
system? First, consider the magnitude of = The location in 
the slab with the greatest deviation in speed relative to Ve 
is either point A or B. In general the orbital speed for a 


ring particle in the Saturnian system is 
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Ry =e 99,790 km 


10,000 km 
— 109,790 km 
Saturn's 


R 
Center | 


mee = 92,000 km =—— iy 7580 km 
m N 
Ro = 104,790 km 


B RING 


Figure 6. A Slab with L = 10,000 km Centered in the 
-- B Ring 
y= Se Ee ee ude 
= oR 


Applying this equation at A, B, and Cc yields 


= -l 
Va = 19.496 kms ; 
-l 
Vp = 18.587 kms ; 
and 
> = 
i O02 5) kins 
AUS , 


all of the particles in the slab orbit Saturn with a 


speed in the range [18.587 mae, 19.496 Ine Ve the greatest 


Bis 


deviation in assigning Mig (equals en! to all of the partveies 
is 2.5%. Now, Consider the @dieectuopee. = The pointven 
greatest deviation in direction, 9, occurs at point E (the 
angle of deviation at point F is equal to that at E)] |p age 
7 gives the computation of 6; the results given there demon- 


strate that the angular error in assiquamd Vo to all of joi. 


~1, 9000 ; 
997120 


6 = 5.006 < 110° - eee 


R, = 99,915 km 


6 = tan 





R /7,| =19.485 kms” 
A v.. =0.975 kms > 
EX 
aie a -l 
eee i Viny =19.460 kms 
m 

Ot 
Saturn 

Figure 7. Calculation of the Maximum Angle of Deviation, 
particles is quite small. The above example verifies the 


CONnCeDE vOL arslanavetleetcy | Vo) for large slabs in theme 
ring (similar results also hold for the A ring). The cCone@egs 
applies even better for smaller slabs. 

Accepting that Ve is a meaningful concept with regard 
to physical slabs in the A and B rings, the requirement that 


it be constant can now be examined. Once defined, the slab 
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feeeleccanslate as shown in Figure $8. Upon reflection, the 


Pree velocity Cannot be strictly constant. Although iy | 


A 
j p 
Slab at 
ame = weak 
iL 
R 
Saturn's A 
R V 
Center 
<_— 
slab at 
time t 
O 
RING 
Figure 8. Translation of the Slab over Time 


equals lv the direction of the velocity vector changes 


pl! 


with time. The angle » which quantifies the change in direc- 


tion is given by 





aes p 
b = ie Gxerere 7°) 
For a slab centered in the B ring (as in Figure 6), 


6 se Sie ose | | (GaeEe: 3) 
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where t 1s in seconds, ¢ in radians. The time required for 


a Variationeonr ——. 745. io 


rad (1°) is 96.1 s. This Gxameam 
suggests that for the A and B rings, the assumption of con- 
stant Ve 1s approximately correct, if the time of passage 

for a molecule through the ring is short. The requirement 
will be met for most passages (i.e., over most y), since the 
thickness of the ring is so small. To see this, then sup- 
pose that: | 


(ly hh -= LO0em, 


(2) R = 105,000 km (molecule passes near the center of 
Ehe Bmime)., 


(3) |¥_| = 20 kms 2,7 

‘=e m 
Over what range of y (recall Figure 2) will ¢ be less than 1°; 
Or phrased differently, over what range of y will the time 


Of passage of the molecule be less than 96 s? In general the 


time for passage for the molecule is 


tc = : (TIL. 
v_,|cos y | 


Substituting the given information, requiring t to be equal 


to 96 S;, and solving for + yretds 


100 


26 5 Sa 
20 < 104 leocua. 


Onhe lowest kinetic energy of a stable orbit (measured 
in a frame of reference fixed with respect to Saturn) for an 
atmospheric particle which originates in the A and B rings, 
as it passes through those rings, is ~20 kms7™ 
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[coca © 10. 


Va ooo, 90.003° . 


Therefore, in this representative case the range of y is 


vo < 160 SOM maa, FOr OO em 1 SOS 


=) 
J A 


virtually any y meets the requirement. 

The requirement that ve be constant, as in the case 
ong vo! cannot be strictly met, since the trajectory of any 
molecular orbit will have some aetna, ot However, the 
extreme thinness of the eine allows the condition to be 
approximately satisfied, so this aspect of the theoretical 
model can be applied to Saturn's rings with no problem. 

That the particles on the slab are fixed with respect 
to Ve Cannot be strictly satisfied, either. Direct measure- 
ments of the relative velocities of the ring particles have 
not been accomplished, but various indirect methods yield 
consistent results: ~l cms + (Weidenschilling et al., 1984). 


Compared to the range of orbital speeds for particles in the 


A and B rings 


Le james < Ses SOs ame : 


(outer A ring) (inner B ring) 


Bote that the gravitational interaction between the 
molecule and ring particles has been neglected, a reasonable | 
assumption since the total mass of the rings is less than 10 
Maat OL Saturn (Null et al., 1981). 


a] 


calculated using equation (II1.6), the relative mouronees 
the particles is insignificant over reasonable passage times 
(less than or equal to ~100 s, as discussed above). Thus, 
in this regard the theoretical model suffers no important 
restriction in applicability Couehnearunce- 
a. oulmacy 

The above discussion has demonstrated that the condi- 
tions upon which the theoretical model is based, are gener- 
ally well-satisfied in the ring ene onnene The only 
potential source of significant concern relates to assumption 
of random particle distribution and the validity of the exponen 
tial .approximation (Section ITI.A.3). Unlike the other asSGmia@ 
tions of the model, these cannot be unamibguously shown to 
hold for Saturn's rings. On the other hand, they cannot be 
shown to fail, and have considerable intuitive appeal. With 
this caveat noted, the theoretical model is adopted for analy- 
Sis of collisions between ring particles and atmospheric 


molecules. 


B. THE UTILITY OF OPTICAL THICKNESS IN APPEDTeAT Tiel 

Having shown that the theoretical model can be viably 
applied to the atmosphere in the vicinity of the A and B 
EINGS> Calculaeronm on Be using the results of Section) Time 
justified. Since particle radius is continuously distriba@eess 
and effectively confined to an interval [0.01 m,5 m] as dis- 
cussed earlier, the collision probability is (from Section 


To) 
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5 
_ 7 -Th 2 
BP. = | eXP Tos eal is 16 siokle hekal amr (Cie 6.) 


: . , — 
when working in a frame of reference with Ws equally eo (ZerG, 


emia is 


-T1h “p “p 2 2 
PL =l ge Pieosmyil — wes") - Ms @ ee) chal aa 


(ET. 8.) 


when working ina frame of reference where Ve is not equal to 
zero. The value for h is well-constrained, while the param- 
eters v_, Vol Gh, YY, ama y° dependyon the location of a mole- 
cule's passage and its trajectory (both of which can be 
determined using Mepel Qeeen mechanics). If n(r) is known 
throughout the rings, then the calculation of P using the 
above equations will be straightforward. Although the func- 
tional form of n(r) is believed to be 

n(r) = nace owe «0-4 
where p may differ over r, its precise expression has only 
been determined for a few locations in the A and C rings 
MiaGoOur et al., 1983). Fortunately, this obstacle to the 
application of equations (II.6) and (II.8) can be overcome 


by using the concept of "optical thickness." 
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Opieme can thickness-“, t, is defined in an operational 
manner: if an electromagnetic beam traverses an lay ome 
particles at an angle uw to the layer normal, and the layer 
has an optical thickness of Tt, then the beam's intensity as 


reduced by a factor of a1 /ESs E 


In general t will depend 
on: particle composition, size, and shape; particle density; 
layer thickness; and wavelength of the incident beam. De- 
tailed analysis of this phenomenon is beyond the scope of this 
paper (see van de Hulst (1981) and Chandrasekar (1960) for 
thorough analyses of light scattering by media consisting of 
small particles). However, the essential attributes of a com- 
plete analysis of 1, can be concisely summaried, generaia 
Following Cuzzi et al. (1969 


For spherical particles ina layer of thickness, h, the 


Optical thickness is given as 
ee) Q_(x,A) mr“hn (x) ar ; (III. la 


where: 
A = wavelength of incident radiation, 


and 


te optical thickness is often referred to as "optical 
depth." The terms are used interchangeably in literature 
concerning planetary rings. Also, since T depends on wave- 
length, one must be careful when comparing the resulesues 
different studies. 
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Ort.) = extinction cross section (a measure 
Geette DOrtion Of incident energy 
removed from the beam by absorption 
eigmamoecattering) « 


When r is much greater than i, then 


imeeeor 1necoemerent incident radiation. 


PP POuTcOnenone Inerdens radlatron . 


Considering visible wavelengths with respect to the predominant 
particle sizes in the A and B rings, i.e., [0.01 m,5 mJ], the 
above values for Q. apply. Therefore, for the A and B 

rings the visible optical thickness when measured with an 


incoherent signal (e.g., stellar occultation) is 
2 
cn au etn (eicie) (alae ll) 


Equations (II.6) and (II.8) can be rewritten as 


2 ae e7t/|cos y'| Grate) 
Ve V 

ne. ee i cos a + (2) -} (TIT.13) 
Cc lcos y| va ra 


Of course, if optical thickness measurements based on coherent, 
visible signals are used to compute T, the resulting expres- 
Sions for Pi WEG atter from eGuactions (LEE.12) and (III.13) 


Ilya factor of 0.5 in the exponent of e. 
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Measurements of 1t for Saturn's rings abound.) !ouaoam 
equations (III.«12) and (IIE.13) to the rings, one musts, 
be cognizant of the importance of the wavelength at which 
tT was determined. Keeping this in mind, analysis of P. can 


proceed. 
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fe coNoee oe ear ols OF THE EXPRESSION FOR 
teewecooorerON PROBABILITY, Po 

fee RESTRICTIONS ON THE PARAMETERS 

When studying Saturn's rings, the frame of reference 
usually adopted features Saturn as a fixed central body 
Orbited by ring particles following (circular) Keplerian 
trajectories. It is also the reference frame of most general 
interest when examining the ring atmosphere, since ballistic 
transport models are based on basic orbital mechanics. For 
this reason only equation (III.13) is analyzed in this 
die 


section. 


Recall that the collision probability 1S given as 


: Vf V 
- - au Bell ya ee 2 
P 1 Sores <p 1 ‘ees a ) 


Mm 


ay (Cig oy) 


Clearly Pes is a function of parameters related to the mole- 
eube's orbit (vir ¥) and those that characterize the 
translating slab ee Note that a and y are not indepen- 


dent, but must satisfy 


cosa + Boca < iar 


3ohat is not to say that the reference frame with slab 
particles at rest is not useful in some practical applications 
(as 1s shown in Section V). Rather, although the collision 
problem may be solved in either frame of reference, the prob- 
lem of calculating molecular trajectories is ultimately ap- 
proached in the frame with v_ not equal to zero, requiring 
conversion to that frame. P 
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Since they are direction cosines. An alternative coordinate 
system which removes this interdependence specifies mole- 
cule's path in terms of y and §, henceforth referred to as 

the "orbital inclination" and “orbital azimuth.” 6 ds "ieee 
defined as the angle between the projection of Vn in the 

x-y plane (the ring plane), and the x-axis (1.e., the directioem 
in which the slab translates). Figure 9 shows the relation- 

ship between § and the direction cosines, the mathematical 


expression of which is given in Appendix D. Employing the 





Figure 9. Definition of the “Orbital Azimueh 


results of Appendix D, the collision probability can be 


rewritten as 


2 


V V 
ee = BS “P . 
Ee =— exPlrass a7 al oie )sin VY COS oA ) } ° (LV 


mM TM 


Equation (IV.1) proves more useful than equation (III.13) 


in the subsequent analysis. 
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Before examining the expression for BP. in detail, physi- 
cal considerations which limit the range of values for the 
parameters are outlined below. Following this, attention can 
then be turned to parameter values of most practical interest 
in the study of the atmosphere of Saturn's rings. 

The optical thickness, tT, of Saturn's rings has been 
studied extensively. Recall from Section III.B that the 
Meesured value of tT is dependent on several variables, most 
feortantly wavelength and coherence of the incident radiation, 
and the nominal size of the ring particles at the location 
where the electromagnetic beam is incident upon the ring plane. 
As shown in Section III.B, measurements based on incoherent, 
short wavelength (e.g., visible) radiation yield values of 
tT that can be most easily used in collision theory (see 
meee dons (III.1l1) and (III.13)). <A stellar occultation 
experiment onboard Voyager 2 provided results that meet the 
above requirements (Esposito et al., 1984). Optical thick- 
ness was found to fluctuate radially, varying from: 

meme 0.) to ~0.4 in the C ring. 

wemee-O./ to ~2.5 in the B ring. 

me ~O.1 to ~0.2 in the Cassini Division. 

eee 0.4 to ~1.0 in the A ring. 
Average tT over the A-C rings is given in Appendix A. Sum- 
Gamezing this information, the nominal optical thickness is: 

mee 0.1 for the C ring and the Cassini Division. 

m) 1.5 for the B ring. 


me 0.5 for the A ring. 
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Subsequently, in further analysis of Saturn's ring atmosonenes 


t will be limited to the inter vas 


The slab speed (which is essentially the ring particle 

speed in the vicinity of the molecule's path through the 
ring) and the molecular speed appear as a ratio, V5/ Ym! in 
the expression for Po: earn ee. the "Speed ratio." Esti- 
mates of H,0 molecule ejection speeds, relative to the source 
particle, for several production processes are: 

62) os Oe! kms + for sublimation (Dennefeld, 1974). 

Cs at | oe for meteroid bombardment (Dennefeld, 1974). 
Estimates for H atom ejection speed are: 


ll 008 kms 1 for meteroid bombardment (Mor£ill et alee 
19s). 


(29. ee 0 kms + for photodissociation of H.O (Durisen, 2o3eae 


2 
Based on the above values, estimates for the speed ratio 
range for molecules newly ejected from particles in the A 


and B rings are 


< 


fOr HO molecules, and 


V 
0.40 < -—F « 1.49 
mM 
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poten atoms, noting that Vin is the sum of ve and the ejection 
welocauty , Nig ies eCeemericity OL Ehe subsequent orbit, 

and the effect of collisions complicate precise specification 
of general limits on the speed ratio, but the above should 
suffice as guidelines. Therefore, further analysis of P. an 


this paper will be limited to a range 


By definition, y and @ are restricted to [0,7] and [-7,7], 
respectively. Considering equation (IV.1), it-can be seen 
femme t -- 

ce = Page) ’ 
and 
Gp Psy) +> Pl(=n-y) = P. (=n +y) 
Cc Ce ome Giz 


Noting these symmetries, no information is lost by limiting 


subsequent analysis of P, to 


and 


= 
eS 
=e 
LAN 

NO] 
= 
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Though Ee will be analyzed over the above intervals, some 
physical restrictions on values for orbital inclination and 
azimuth are discussed in the following paragraphs. 

When a given production process is considered, the initial 
limits on y can easily be determined as shown in Figure 10. 


For example, the y-range for H,O is: 


ay, ye , 
: \ 
a 
ny . ay 
Pp 


\ rd 


/ 
4 


LOCUS Ob OSs tele 
ejection direction 


(a) (b) 





Figure 10. (a) General Relationshw) bee —a venk V 
and Vy. 
(b) Case eae Minimum Value of y 


(lie “89%. 4s) ar = 293 for sublimaruone 

(24 “286360 Yess. cy: Boe oo for meteroid impact, 
Similarly, for H atoms, y is initially contanecciiae. 

(1) -79.6°2.-= “yy “= S00 for meteroid impacr: 


(2)' 5320° 5 46" S452 7 aoe for photodissociation of H,0O 
on a particle surnaeee 
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As collisions occur between molecules and ring particles, 

the range on y will change. If a collision is assumed to be 

inelastic, then upon impact the molecule will rebound with 

reduced kinetic energy relative to the particle. The overall 

effect is that Y should go to =m asemone  collascions Occur. 
For atmospheric molecules originating on a ring particle, 


the initial orbital azimuth will certainly be "forward- 


erected," 1.e€., 


MN} Fr 


as long as the ejection speed is less than the particle speed. 
This again follows from the fact that Va is the vector sum of 
v, and Voss Subsequent collisions should not alter this 


meecertion. 


ID. BUNCLEONAL DEPENDENCE ON OPTICAL THICKNESS, Tt 
Pa monotonically increases with increasing tT, as can be 
clearly observed in equation (IV.1). The rate of increase 


memciven by (using equation (IV.1)) 


dP V V V 
=f : *D) _P)¢ 2 (2) s:; 54 (2 
Bese alt a ) sin 7 COS 6 + ( o ) “expt ssey7 = ) sinycos8+ Ce ) 


m m mM 


ae 


Here the effect of the other parameters on the rate with 
which a exponentially approaches one is explicitly displayed. 


Appendix E contains examples of the collision probability's 
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dependence on Tt given selected values for the speed ratio, 
orbital ‘inclination, and oreittal azimeeu 

The physical interpretation of this functional behavior 
is simple. Recall that for measurements at appropriate 


wavelengths, 
2 
T = f pagel (se) ole se 


If h is held constant, then higher values of tT imply greater 
number density, n(r). Intuitively, one would expect the 
collision probability to increase with increasing number 
density (all other things held constant), which is substan- 
tiated by the mathematical expression for Poe Similarly 

1f n(r) is held constant, increasing t implies that h mise 
increase, the length of a molecule's path must incre 


thus increasing its exposure to collision. 


V 
C. FUNCTIONAL DEPENDENCE ON THE SPEED aire —~ 
m 

The behavior of the collision probability as a function 


of speed ratio is very interesting. Figure ll gives a repre- 
sentative example of the variation of P, asa Functionves 
ee (see Appendix E for the results of more example calcu- 
lations. For 6 less than sm, a minimum for ae OCCUrS ae 
V/V greater than zero. Why does this happen? 

Some insight into this phenomenon Can be gaincame 


considering the partial derivative of Pe with respect to 


Val Ts i.e., 
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Figure ll. Collision Probability as a Function of the 
Speed Ratio for tT = 1.0 and y = i/4n and 
Various Values of 6 
se Y “p. —T ue Yo 2 
= = Teos yp 2sinvoosé+2 (7) lexpl tess 7 a a ala 
a (<E) 
m 


ive) 


ey 


With all variables except Ven) in held ceénstane- Pa has a 


minimum at 


= “SIN yeecos oe (IV oge 


abo 


From Appendix D, the right hand side of equation (IV.3) is 
identically equal to cos a, so expressing that in terms of 


the direction cosine gives 


Vv = ee Cosmo (IV.4) 


a minimum of P. occurs when the component of Vn in the direc- 
Cero a 1s equal in magnitude to vl. Notice that relative 
tO Variation in eae Pp. Cannot assume a minimum value for 
speed ratio values greater than one. 

Since the result of the above mathematical analysis may . 
be difficult to reconcile with intuition, an example wis 
given to clarify its physical meaning. Suppose that a blind 
mouse lives ina hole in the middle of a single-lane, one- 
way road. Randomly spaced steamrollers traverse this road, 
all moving with speed VS The mouse, smelling cheese in 
another hole a few meters down the road, dashes for that 
hole with speed Ves (in the same direction as the steamrollers). 
Given ve fon whack ws will the mouse have the best chance 
to reach the cheese unharmed, neither being overtaken by a 


steam roller, nor running into one? The answer is that the 


a2 


PeObantt ity of collision is least (in this case zero) when 
Vo equals Vi Bieta a greater than Wier steamrollers will 
tend to overtake the mouse; as Ns gets very large the proba- 
eelity Of collision goes to one. For Wis less than v1 the 
mouse will tend to run into a steamroller. The mouse's 
situation is analogous to that of a molecule passing through 
mearng. che probability of collision for the molecule is 
minimized when it "goes with the flow" of the particles as 
meen 25 1S possible. This minimum value of an Wie not in 
general be equal to zero for molecules passing through the 
ring, as it is for the mouse. This is the case, because 

the molecule crosses the ring plane of some angle (Speci- 
med by y).- 


A final comment on the influence of el he on Po is that 


P > ae 
S 
as 
V 
wae > oa 
vm 
regardless of the values of the other variables. For ve 


very large, the ring loses its discrete character becoming 
essentially a continuum; collision becomes inevitable. 
Pee FUNCTIONAL DEPENDENCE ON ORBITAL INCLINATION, y, AND 
ORBITAL AZIMUTH, @ 
y and 8 specify the path of a molecule through the 
ring. Figures 12 and 13 show two examples of the 


oe) 
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3 = =r 
5S Me ) 
0.9 9 7 oul Note: Ea 1s 
ul undefined at 
= i, 6 = 0 
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l; 
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Figure 12. P as a Function of y and 6 for tT = 1, 2 2 
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0.9 
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Figure wise Ee as a Function of y and @.fonueee-e- = = 7 
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functional dependence of Ee Sieeroweal Inclination and 
eaimuach . 
Bicse, consider Geaesbehavior Of Pa with respect to 4. 


The partial derivative with respect to this variable is 


oP 
=a = TTP) sinvsin 1-2 <E sinyooss 
ae -P) 53 “py? 
“ om ] explsszy 1-2 ca ina } 


>) 


From equation (IV.5) it can be deduced that Ee 1s minimized 
when--6 equals zero and maximized when 8 equals 1 (when all 
other variables are held constant). As expected, both Figures 
[i2eand 13 clearly confirm this result. Also, notice how 
rapidly se approaches its maximum value with increasing 6. 
The story of the blind mouse from Section IV.C provides the 
physical interpretation for this behavior: all other param- 
eters held constant, collision probability is least when 
going with the flow of particles (or steamrollers), and is 
greatest when moving against them. 

The variation of Bx with y is a little more complicated. 


The partial derivative of P with respect to y is 


oP V v5 9 eZ 
= = [t secy tany(1-2()sin ycos 6+(*)*) 
Y Vv, be 
m m 
v5 Vp V5 5 -l1/2 a v5 V5 > 
- t(—)cos6 (1-2 (—-) sinycosH— { eet Sy COST Ha) | 
om a Y ae lexptycos y] wa u a 
(1 Veo} 
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From equation (IV.6), the condition for minimum Po, when 


all variables except y are held constant, is 


VV 
(—P) cos 6 
Sec 7 cata aa m 
2 V : (iV . 7a 
1 +2tan y * (PB) ¢ 
vin 


This imposing expression has a simple physical meaning: 
(l) dif v /Vn is less than one, P. is minimized whe 
the_flagnitude of the component of vp, in the direcemar 
of v_ is equal to v 
p p 
(2) if vyo/Vvpm is greater than one, Pe 1S minimized when 
the magnitude of the componen ager V5 in the direction 
Of vp is equal to vp- 
Once again, the blind mouse analogy may be helpful in under- 
standing this behavior. 
Before leaving this section, a few words are in order 
to clarify possible concern over the point of "undefined" 
collision probability shown in Figure 12. This anomaly venus 
occurs when Ney iy is equal to one. For this particular speed 
ratio, when 6 is zero, and y is 7/2, then the molecule is 
moving in concert with the ring particles. Thus, 22 eee 
not initially in contact ewmeaea Darencler it never will be. 
On the other hand, if initially in contact with a paneveiee 
1t will never separate from that particle. Therefore, the 


notion of a collision probability is meaningless under this 


unique circumstance. 


Vee oO LSOCUSSTON 


A. GENERAL COMMENTS 

To effectively exploit the collision theory developed in 
this paper, one generally must couple it with a ballistic 
transport model which provides values for the speed ratio, 
orbital inclination, and orbital azimuth. More accurate 
eee mation of atmospheric density will result, although 
improvements in this area are limited by the lack of knowledge 
of details of a molecule-particle collision (i.e., the stick- 
maegecoefficient described in Section I). The frequency of 
collisions also affects the spatial extent of the ring 
atmosphere. The trajectory of a molecule is initially 
determined by the nature of the production process, while its 
Solution is dictated by collisions. 

Consideration of the special case of low energy, iso- 
tropic production processes (e.g., sublimation) allows one 
to apply the collision theory in making inferences on the 
Spatial extent of the ring atmosphere. 
feo ePLICATION TO AN ISOTROPIC, LOW ENERGY PRODUCTION 

EROCESS 

Define an isotropic, low energy production process to 
meeeemechanism which ejects molecules from ring particles 
aeem that; | 

(1) the direction of ejection relative to the particle 


is random. 


De 


(29) See <a Woe which implies Sol Gi —— 
The orbit of a molecule, prior to collisiGns 772i 
particles, will be very nearly circular with a radius eames 
to that of its parent ring particle. The "Gretta 
will be only slightly inclined with respect to the ring 
plane, e.g., for sublimation the orbital plane is incline 
at most 1° (see Section IV.A). 

When a collision does occur, two outcomes are possible: 
the molecule "sticks," thus removing it from the atmosphere; 
or the molecule rebounds, assuming a new trajectory. In 
the latter case, what can be said about the new orbit? 
Since. the collision will surely be inelastic to some degree, 
the speed of the rebounding particle relative mG the ring 
particle is less than its incident speed. The speed ratio 
remains almost equal to one, implying a circular orbit. 
Following the approach used in Section IV.A, especially 
Figure 10, the primary effect of the inelastic collision is 
to reduce the inclination of the orbital plane. Thus, 
collisions act to reduce the spatial extent of the ring 
atmosphere; the details of this reduction being determined 
by the frequency of collisions, and the actual inelasticity 
Of an, andividual collision. 

The more frequent collisions are, the more closely cone 
fined to the vicinity of the ring will be the ring atmosphere. 
The collision theory can be used to examine collision fre- 


quency. Specifically, the average probability, as a function 


26 


Peete, tnac a molecule has not suffered at least one 
collision is calculated below. 

When a molecule is initially ejected, it may collide 
with a ring particle before ever leaving the ring plane. 
mo, calculate the probability that it collides with at least 
one ring particle, adopt a frame of reference in which 
the particles are fixed. The probability that the molecule 


is ejected from a particle in dz' is (see Figure 14) 
Ci ee ec) ce, (yea) 


while the probability that it is ejected into solid angle 












aQ is 
— it 7 J ! t 
P(Q)dQ = x Siny Cow Sa Si (Vi223 
yp 
ip 
| 
"Effective 
BOLE h NESS | Ejection Angle 
Figure 14. Ejection of a Molecule from a Particle at z' 


32, 


Since the production process is isotropic. ~The probae sae 
density for a molecule ejected at y' from z' colliding with ‘ 


a ring particle before leaving the ring 


= =( (he2ee ' 7 
1 expip Seen Ca es OS 
q ' = 
+T-Zz' TT 7 
1 XRF cosy oe 


which follows directly from equation (II1:.12),) whenmege 

effective slab thickness is taken into account. The aver- 
age probability. of at least one collision for a molecule 

produced by an isotropic process 1s found by integrating over 


the product of equations (V.1), (V.2) anaes aoe 


1 
Z h 
= 7 1 : eco -t(h-z') ae 
Re J iF x hsin y' [l expt Pcogyt az dy 
ei 
a Se Beenie zy ; 
+ Ms - a hey a Lt explicogy 142 dy 
pu 
Making a change of variables, x = cos y' and xX = - cosy’ 


in the first and second integrals, respectively, then inte- 


grating, yields 


(V 245 


| 
I 
= 
| 
al ae 
— 
© — 
x 
= 
i 
49) 
h 
4 
se 
* 
QO, 
x 
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Therefore, the average probability that a molecule initially 


(i.e., t = 0) escapes the ring without collision is 
Pp = 1-Pp = Li x(1 - e /*) ax (V.5) 
esc & T “9 ‘ 
able 2 gives Bae Pome eiviec@or et On emiiterest for the rings’. 


Notice that there is a significant chance of collision immedi- 
ately after a molecule is ejected; suggesting that the spatial 
extent of the atmosphere is biased toward being closely 
confined to the ring even before subsequent evolution of the 


system is considered. 


TABLE 2 


AVERAGE PROBABILITY THAT A MOLECULE INITIALLY 
(t = 0) ESCAPES THE RING PLANE WITHOUT 
SOVLEDING WEEH RENG PARTICLES 


P 
T esc 





During each orbital period, a molecule in a circular 
Orbit traverses the ring two times. If that molecule has 


Meme cOllided with a ring particle, y‘ retains its initial 
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value for each subsequent passage of the rings, and the 
probability that it does not collide with a ring spar cee 

: . 4 -2T > hal : ‘ 
during each period is Ss Gree TT. The probability tia 


the molecule is ejected into solid angle di without coliaaaae 


with a ring particle is: 








ee <7" = 


ee cs -eplisen Z 


i8og T; is the number of completed orbital periods, then the 
probability that the molecule has not collided with any 
particle as a function of time, given in terms of 75! 1s 

gosy” 7 i 

(i: explaoat ——7}, 0<y' <n 


TT ~2tT,, 


oe 
esc'Th) = J. 3siny exp (Tess yo 


ae 1 
=e oe Piet ism <y' <n 
Owe 
alt =—2 /X 7 
P(t) = =f xe (santo (v.6) 
esc’ p a, 
14 


This result can be easily understood by considering 

a slab of thickness 2h. Extending this PeéasOmuggp ee 
periods have been completed, the probability of no colliszea 
for a molecule initially escaping is exp{-2TT »/|cosy' lee. 
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Table 3 gives the probability that a molecule has remained 
free of collisions,as a function of time (given in orbital 


periods) for various optical thicknesses. These results 


TABLE 3 


ans : 
VARIATION OF Baas ie te AND T 5 


esc 
"s i = (Os 1 Ose S (cue — oe! ae 6. tT = 1.5 =a 
0 0.837 a BBE 0.390 0.296 0.235 
1 a. 51e 0.106 0.021 0.005 0.001 
2 0.357. | 0.028 0.002 iL. Besar 4 ~0 
3 0.255 0.008 1 @asaiie Te 26 0 
4 0.186 0.002 =o 0 ~0 
5 0.138 0.001 ~0 20 ~0 


Seeongly imply collisions are quite frequent for isotropic, 
low energy production pees oseecr TiS ein eeurneaimplies that 
the ring atmosphere associated with this process will tend 
to be more concentrated in the vicinity of the ring than 
Smalysis Of the distribution of particles prior to collision 


would indicate. 


faoN MULTIPLE COLLISIONS 
The collision probability theory developed in this thesis 


Cannot be used to precisely calculate the probability of 
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multiple collisions, i.e., the chance that a molecule col- 
lides with two, three, etc., particles during a single 
passage through the rings. Computer-simulation models may 
provide the best way to approach this problem. 

Based on the results of the preceding section, the 
probability of at least one collision is guite high in the 
A and B rings. This suggests that multiple collisions eau 


be a significant factor in analysis of the ring atmosphere. 
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Vin ~ CONCEUSIONS 


In this paper, a model for collisions between ring 
particles and molecules in the ring atmosphere has been 
developed. The resulting expressions for the probability 
of at least one collision when a molecule passes through 


the ring are 


t's 
P. = l- 22 ress sup ’ eet. 2) 


for a frame of reference where the ring particles are fixed, 


and 
ee re 2 
PL, = l-exptyess-y7 nee Samy eel ie (ye) 


for a frame of reference fixed on Saturn in which the parti- 
cles are seen to orbit the planet. They are more complete 
than those currently employed in analysis of the ring atmos- 
phere, Since specifics on molecular trajectory and the rela- 
tive velocity of a molecule with respect to the ring parti- 
cles are included. Thus, this analytic collision theory is 
compatible with computer-simulated ballistic transport 
models, since they include similar effects. 

The frequency of collisions influences both the total 


population of molecules in the ring atmosphere, and their 
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spatial distribution, as shown by Dennefela (1974 cae 
knowledge of molecular sources is still evolving, and uneew. 
standing of the dynamics of an individual collision (eegee 
Dennefelid's “sticking coefficient") is in Pesan eames 
improved estimates of molecular populations cannot be made 
based on this revision of collision theory, alone. Infer- 
ences on spatial distribution can, however, be made. When 
a low energy, isotropic production process is considered for 
optical thicknesses characteristic of the A and B rings of 
Saturn, collisions are found to be guite frequent. This 
implies an atmosphere closely confined to the vicinity of the 
rings. Extending this result to more energetic production 
processes (e.g., see Section III.A) implies that the ring 
atmosphere will be toroidal in shape. Thus, assuming a 
toroidal shape will probably be more productive than assum- 
ing ee eiCal shape in future studies of the ring atmosphere. 
As a final note, the collision theory used here in con- 
nection with the ring atmosphere should be useful in another 
related application. Ballistic transport of macroscopere 
particles produced during meteroid impacts may account for 
some of the radial structure in Saturn's rings (Ip, 1983; 
Durisen, 1984). These particles form a "chip halo" in the 
vicinity of the rings. As long as the chips are small vega. 
tive to ring particles, the model developed here applies 


for calculation of collision prebabige emcee 
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APPENDIX A 


RING NOMENCLATURE AND DIMENSIONS 


Saturn's ring system consists of seven distinct rings. 
The A, B, C, D, and F rings are characterized by distinct 
boundaries. Their nominal thicknesses are not more than 
=—20'0 m (Cured et al., 1984), and may well be less than 5-10 
m (Bridges et al., 1984). The E ane! © rings are diffuse 
structures with poorly defined boundaries and nominal thick- 
nesses of 210" and Las km, respectively (Cuzzi et al., 


1984). Table Al is a summary of the radial structure of the 


ring system. Note that the F and G rings are quite narrow. 


TABLE Al 


RADIAL STRUCTURE OF SATURN'S RINGS 


Boundaries 
Ring/Region Ro km mass 
D SLs sile 2o6 Bey 70274) 510 ? 
C 193551595 Fd Sieeee woe 2x10"? Ms 
B 1.525-1.949 92,000-117,580 5 x10° Ms 
— io = 2 400.5 LLG, BOA Ske 1x10? Ms 
A DROPS Mw 122,170-186,780 110° Ms 
F 7. OE 140,180 2 
G 2.82 170 100 0) Eeioe! Ms 
E 3-8 181,000-438,000 ? 
Re —wo0,330) km (Saturn's radius) Me = 5 GOB sao. kg sa il 


ow 


The Cassini Division is a well-defined gap between the A and 
B rings, which 1s by no means | emer es 

The optical depth, T, is a dimensionless parameter 
directly related to mass density. Table A2 gives average 
optical thickness for the main rings. Optical thickness in 
the other rings is much lower. Since optical thickness 
increases with increasing mass density, it follows that the 
A and B rings should be the most significant sources of a 


ring atmosphere. 


TABLE A2 


AVERAGE OPTICAL THICKNESS OF THE MAIN RINGS BASED 
| ON STELLAR OCCULTATION OBSERVATION GS: 2G 
VOYAGER 2 (ESPOSITO ED Ai ee 


REGION BOUNDARIES (R_) “AVE 
inner C 1.24-1.39 OF0'S 
Gucer 7 SS ele 52 Oi 
inner B 1.52-1.66 Lik 
middle B 1.66-1.72 linis6 
outer B ipo. 9.5 Le 
Cassini Division 9 S202 Ore 
inner A 2202226 Oc 70 
outer A 2 Oe 0 238) 
Re = 60,330 km 
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ae Nea Ss 


DERIVATION Oper. WHEN Vo IS NOT EQUAL TO ZERO 
i SS eee 


In Section II of the main body of this paper, an expres- 
sion for the probability that a molecule passing through a 
ring collides with at least one ring particle, aa was 
derived in a frame of reference in which the ring particles 
are fixed (l.e., “ equal to zero). Here it will be demon- 
strated that the same expression for PB. results when the 
derivation is made ina frame of reference where the slab 
translates with v5 Netwequals ce Zero. 

The model upon which this derivation is based is the 
same as that described in Section II. Reviewing, a molecule 


passes through a region in space which contains randomly 


distributed (spatially) spherical particles. This region, 
the "slab," is of large extent, L, in two dimensions- and of 
relatively small thickness, h. The entire slab moves through 


PY) are r) ry ° 
space with constant velocity, Vi! timo -eectlOn nOrmaL EO 
the slab's thickness (V7, is measured relative to some iner- 
tial frame of reference). The radius, or "size," of each 


particle is restricted to an interval [rs 


eee! Wich or 


min 
much greater than the nominal radius of the molecule, and 

with a much smaller than h. The size distribution of the 
particles may be either discrete or continuous. The number 


of particles in the slab is further assumed to be large. 
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Finally, the number density as a function of size is denoted 


by n(r), defined as the number of particles: 


(l) of radius r per unit volume (discrete size 
distributions). 
(2) ain the radius interval [r,r+dr] per unit volume 


(continuous size distributions). 

What is the probability that the molecule collides with 
a particle in the slab? In Section 11 this )quese oie. 
answered with relative ease by converting to a frame of 
reference in which the particles are at rest. Below, no 
such simplification is made. 

When Vo is not equal to zero, the notion of cylindrieas 
"colfFision regions," used an Section II, loses its validity. 
Another approach must be taken. Suppose that the slab is 
subdivided into regions of very small thickness, Ah, and . 


finite area ee (see Figure Bl). If the probability that the 


baa 


_——$—$—— ee 


Figure Bl. Edgeview of the Slab and a Representative 
Subregion 
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molecule does not collide with particles of size r whose 
centers lie in a specific sublayver can be calculated (call 
this Pig fe ARs), Piewmene total ~robability of no colli- 


sion with particles of size r will be 


Since all of the small sublayers are mutually independent.*? 


mie Ccollasion probability, following Section II, is 


for discretely distributed r. 


The problem now is to find an expression for F ae age 


First, the case of discrete r is solved from which the con- 
tinuous solution can be deduced. Adopt a point of view in 
which the molecule has radius r, and the particles are point 
masses (on reflection this is analogous to the real situation 
as far as collision probabilities are concerned). Then, 


Pig (F471 4A) can clearly be seen to be equivalent to the 


probability that the volume in the subslab "swept out" by 


two observations on equation (B.1) are in order. First, 
there are infinitely many factors in the product over i. 
Second, Era ao rena) must be the same for all i, since in this 
model the statistical properties of each subregion are iden- 
mea Equation (8.1) may be difficult to interpret at 
this point, but its meaning will be clarified later in this 
appendix. 


qk 


the passing molecule has no particles centered within it. 
The probability of no collisions with partie ves sense 


Yr. in subslaben] as. 
5 1 


L“Ah.-(Vol. Swept Out) n(r,)L7sh. 
Pio (t5rAhy) = (~~ =] - =? (B.3) 
Lah, 

Since the particles are assumed to be randomly distributed 
throughout the slab (see Section II for a more detailed 
explanation of this assertion). In general, the mathematics 
involved in finding this "volume swept out" is extremely 
tedious. To demonstrate the method, the special case of y 
aaa 8 zero or 7 (molecular path normal to the plane of 
the slab) is solved below. 

On what portion of its path can a molecule potentially 
collide with a ring particle of size Bs al) 9 Ah .? Consideration 
of Figure B2 shows that the pertinent part of the path is of 
length ak If time, t, is set equal to zero when the mole- 
cule initially encounters the subslab (Figure B2(a)), then 


at any t such that 


the molecule will collide with any particle is a disk of 


radius (call this the “inte@ractrven esse auc 


_ Shae ee 
R(t) = (2rv_t vat ) ; (B.4) 


qZ 


Molecular 
Path 


eh. 
1. 





Figure B2. Relationship Between the Molecule and a 
Subslab When Collision with Ring Particles 
Centered in the Subslab is: (a) First 
Possible (t = 0); (b) At Any Time, t, 
Suen that O-< £ < BE ad 


where RL is defined as shown in Figure B2(b). As the inter- 
action radius evolves, the particles in the subslab move 
with speed eos Considering these two factors together, it 
can be deduced that an oblong “interaction region" in the 
subslab must be particle-free, if no collisions are to occur. 
To actually calculate this volume, first relate Ry to the 


distance that the slab moves in time t, call this Dei clearly 


Dy ic) nn Oe (B.5) 
p 
Substituting from equation (B.5) into equation (B.4) yields 


Re) ee) (pa) (B.6) 
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Figure B3 is a sketch of equation (B.6). "The eblengeeness 


in Figure B3 (bounded by the graph for R_ and its refle@umes 


JE 


across the D.-axis, which is shown by a dashed line) cannot 


Ry (Dy) 





Figure B3. Rt as a Function of DS 


represent the projection of the interaction region into the 
plane of the subslab. Consideration of instantaneous inter- 
action region, represented by the circle in Figure B3, illus- 
trates that the actual interaction region is larger. Knowing 
R, (Ds), the function enveloping the interaction region can 


be found. Call this  functacn Ro. 


Consider Figure B4; it shows one-quarter of the inter- 
action volume's projection into the plane of Ah, . At a given 


Gime <t) 7 De. equals Vint |: and the interaction radius is RL. 


The value of R" (D>) for any D, is 


| 2 rae VON ils so ivan 
({R, (D7) l= |D.. pe | <Roe 
Rt (Dt) = (B.7) 


| 0 - otherwise 
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R>/RE 
Gee ee —_—— ' 
oa | 
” 
oo 
R 
= - Jt 
a 
y f | 
Se aa 4 | 
a, ‘ | 
Pa ‘s ' 
7 ee esc tp: | 
7 ‘* i Bey 
4 ! NI 
D 
. ' V ale 
Pr (min) °F x oe ey 
oT t - V 
Ry {D,) J Ym 
Figure B4. Geometry for Calculating RI 
RI (D,) is the maximum value of RU (D>). Setting the first 


derivative of Rt Wien 6esoce eaeSo D. equal to zero, and 


peeving for D. gives 


i = Pp 
D = ee a ey a -f (B.8) 
2 


meen equation (B.6) is substituted into equation (B./). 


S@ecettuting from equation (B.8) into equation (B.7/) gives 


the maximum value of R,(D,) for a given Des 1.e., 
V 
Z Zoe ant 
, vo es cam 1/2 
R,(D_) = ae) es a) : (B.9) 
P 1+ (—) 
Pp 


Using equation (B.9), the region in the sublayer Unaemuee 
be devoid of particles of, size e if no collisi¢6cn tees 


occur between the molecule and such particles, has 


2) 
1g a 
Volume = (Ah,) 4f R*(Dj) aD. , (B.10) 
I (min) 
where 
Vv ily 
at Dye = _p)? | 
meee [ (—) {l + (es ae ese : (B. le 
Mm mM 


When the integral in equation (B.10) is evaluated, 


NGA 


a) AR (Ba Pag 
a: 


2 V 
Volume = mr tl i (—E) 
m 


and substituting into equation (B.3) gives 


V Ae A 
ee ree lee a2) Se Ane nee ae 
ae 7 ie a J 1 
ee = | Se : (B.13a) 


bw, 
a 


Equation (B.13) can be rewritten as 


V yy 
nron(x,) (1 — 
: at (B.13b) 


P 9 (F./ Ah. ) =~ oe 
J 2 L ele) 


Nn 
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Summing over i yields the probability that no collision 


occurs in the entire salb, l.e., 


AZ 
2 Vo 2 
mae SG am } een 
feeir.) = {1 - ] (B.14a) 
no 2 
emery) 
J 
Multiplying the second term in the brackets in equation 
(B.14). by h/h yields 
dyad 
2 “ep. 2 
oe ee 7 om? teat 
eer.) = {l1- hi aa J (B.14b) 
J L“n(r.)h 
J 
but the total number of particles of size Bo 3 n(r,)L°h, 1s 


assumed to be large, so equation (B.14b) can be rewritten 


in exponential form as 


= nee “P 
P (r5) = exp{ moan(r5){l ea) } h} (B.14c) 


mM 


fitjererore, the collision probability is (from equation (B.2)) 


D V9 ee 
P. = l- I] expt-mrn(r.) {1 + (=) } ela 
J m 
D ies) 2 
ae = l1l- exp{-1() ran(r,)) {1 + (>) } eee (Besa } 


J) mM 


or - discrete, or by extension 


ia 


> aie, Vy 2 
P, = 1 - exp{-n(f r°n(x)dr){) @= (B.15b) 
mM 


fOr CONEInuCuUS sa. 


Comparison of equations (B.15a) and (B.15b) with eCquamaams 
(II.7) and (11.8) (with y equal to zero), reveals the equiva- 
lence of solving the collision probability problem in two 
different frames-of-reference. That the results are the same 
1s not surprising. What is interesting is the relative 


difficulty of the two approaches. 
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APPENDIX C 


PeoOME MOmen err RENEE CONVERSION 


iimoceection Ii Of@this thesis, the collision probability, 
Ee? is calculated in a frame of reference in which the ring 
particle velocity, vos is equal to zero. When coupling the 
collision theory with ballistic transport models, a more 
convenient frame of reference is a localized cartesian sys- 
tem fixed in the vicinity of the intersection of the mole- 
cule's path and the ring. The x-axis is in the direction of 
Tonnor the ring particles; the y-axis is directed radially 


Outward from Saturn; and the z-axis is normal to the ring 


plance (see Figure Cl). Here the direction cosines of a, §&, 





Figure Cl. A Coordinate System Fixed at a Specific 
Location in the Rings 


yee, 


and y quite accurately describe the path of the molecule 
through the ring which is essentially linear owing to the 
great radius of curvature of the molecular trajectory 


(~107 km) relative to the extreme thinness (~10 71 


km) of the 
Fangs. 

Let the frame of reference, for which the ring )pan eee 
in the region through which the molecule passes are at rest, 
be called the "primed" system. The simplest coordinate 
system that can be used to describe the molecule's path in 
this frame of reference will have the x'-, y'-, and z'-axes 
in the same direction as the corresponding axes in the un- 
primed system with the y'-axis in the direction of va (see 


Figure C2). a', 8', and y' yield direction cosines deéScuiiume 


the molecules path in the primed system. 


Z 





Figure C2. Direction Cosines in the Two Coordinate Systems 


In Section II the collision problem is attacked using 


the primed system. The expressions for PE derived there 
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(eqeei@nvcmGhl. os) cand (fi-6)) Contain a factor of |cos a. 
Piamexprhession £Or Ehas factor in terms of a, 8, y, ee and 


v_ must be found to represent P in the unprimed system. 


First, note that: 


> > 
i rt ee (ele 


gives the relationship between the molecular velocities in 
the two systems. Writing vector equation (C.1) in component 


form yields: 


ae = Vv —-— Vv ;} (Er Za) 
mx mx p 

Va = Vv : (C225) 
my my 

i Ao | (eee) 


keeping in mind the definitions of the two systems given in 
Megane Cl. Equations (C.2) can be rewritten in terms of 


moe aQirection cosines: 


CO sme rr 9 Vou COS .a@. = Vv. 3 (G2 3a) 
m m D 

t ' — 
ve cos 8 = es Gos 6; (Cie)) 
es COs Se=— Va COS ¥ (Gea) 


eal: 


Squaring equations (C.3) and adding yields: 


vi* (cos“a' Poosak Peose ae = v- (cos‘a Beosee Osan 


+ v - 2V Vv COsm (CC. 45 
Pp m p 


Since the sum of the squares of direction cosines is iden- 


tically equal to one, equation (C.4) becomes 
2 
V SD Me Ma i SO Sie. (Cam 
P m Pp 


Prom equavlen A(ee3c) 


i, os cos y 
Vin oie errs ve S (Com 
which, when substituted into equation (C.5), yields 
2 Gee me = vie + a — 2} .ViesCosnomse (Cr 7 
m°'cos y m p tie 
| yes 
Solving for |cos y'| gives 


Meas 


1 z l nee “Pp 


Equation (C.8) is the fundamental result of this derivation. 
When substituted into equations (II.5) and (II.6), the result= 


ing expressions for ee equations (II.7) and (II.8), have 


See 


the advantage that they are written in terms of the frame 


of reference most often utilized when studying the rings. 
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APPEND Ia 


RELATING ORBITAL AZIMUTH TO THE DIRECTIONS CCl ea. 


Consider a unit vector, u, as shown in Figure Dl. Its 
direction can be specified by the direction cosines (a,8,y), 
or by (y,8) where 6 is called the “orbital azimuth." How 


are these two representations related? 





“W 


Figure Dl. Angles Used to Specify the Direction of u 


First, note that 


eoeee + eee + cosa = ies 


SIG 
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nN 


The x and y components of u are 


u — eos Ooo; CD 2a) 
x 
i =—mnGc@s 'ba.. (Dee2ion 
v7 
By definition 
u 
tan@ = —, Coss) 
a 
x 


and substituting from equations (D.2) gives 


cos 8 (D.4) 


Ean G 
cOoS a 





Squaring equation (D.4) and substituting from equation (D.1) 


yields 


; _ 2 ae 2 
tan G = EE . (Da) 


cos a 
Rearranging equation (D.5) gives 
(Ll + came Coe aa a le cosy ; 


put 


1 + Semace = aae-k and 1 - 60s-y = etna ; 
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so 


2 2 ; 
sec 90 COS @ 7 = secu s, 


Finally, 


COS G@  =""sin 7 scoceoa 





which is useful for expressing P= in terms (ote 220s 
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Pee Slibin = 


Se een ocleroeor CALCULATION OF P 
—S Oe ee 


This appendix consists of several examples of results 


fetmetie Calculation of ae based on the equation 


V V 2 


= zs ail = 9) SE) 6; we 
. IL expl7aos yT IL ae) 22 r.eOs 8 + (a ies 


mM 


Calculations were made over the range of optical thickness 


Seeeinent to Saturn's A and B rings: 


The figures have been arranged to allow one to compare the 


effect of varying the parameters of the equation. 
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Faqure, tebe 


Collision Probabiviagty Sto cee, 


8 8 


Figure E2. 


CoEyoronerPrOooabilaty for + 
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APPENBI 


ANNOTATED BIBLIOGRAPHY 


The purpose of this appendix is to briefly review seme 


papers which should be of particular interest to one who 


wishes to study ring atmospheres. 


General Theory 


a 


Blamont, J. 1974. "The 'Atmosphere' of the Rings of 
Saturn," in NASA SP-343, edited by Palluconi and 
Pettengill, 125-129. 


Blamont's paper is the transcript of his comments, 
delivered at a planning conference for the Pioneer II 
mission. His theoretical model is not detailed but 
results for the H70, H, and OH atmospheres are ug im 
Estimates of Lyman-a emissions are also presented. 


Dennefeld, M. 1974. "Theoretical Studies Of an )@iiaee 
phere around Saturn's Rings,” in Exploration of the 


Planetary System, edited by Woszczyk and Iwaniszewska, 
471-481. 


Dennefeld was Blamont's student. His paper predicts 

the density of the ring atmosphere (and that associated 
with Titan). The sources considered are sublimation, 
bombardment by meteoroids, solar and interstellar wind 
effects, and photodisssociation Of H>@son eEhe ime 
particles. The resulting estimates are outdated for the 
most part, primarily because the production mechanisms 
are now bettertmderstood. The conceptual model for 
calculation of atmospheric densities, however, is of 
enduring significance. It concisely gives a qualitative 
way to describe the atmosphere. One questionable asser- 
tion made in the paper is that the H atmosphere will 

be spherical due to collisions between the H atoms and 
ring particles; this is at odds with the rmesulesme. 

the present thesis. Also, Dennefeld estimates the 
collision probability to be equal to theveperc as 
thickness. 
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Ballistic Transport Models 


1. Durisen, R.H. 1984. Plsanmspone ELrects due to Particle 
Erosion Mechanisms," in Planetary Rings, edited by 
Greenberg and Brahic, 416-446. University of Arizona 
Press. 


Durisen presents a detailed model aimed primarily at 
Peccilceiormeom  Ghip Haver’ 1.€.,"Macroscopic ejecta 
resulting from meteroid impact, formation (the model is 
in general applicable to the ring atmosphere, too). 

The assumption of “prompt absorption” (1.e., particle- 
molecule collisions are inelastic) is made for calcu- 
Tae tonse tai nis paper aesin calculating collision 
Buooaoutty, Durisen accounts for “slant path” of a 
molecule through the rings, but it is not clear whether 
he considered the relative speed of the particles versus 
the molecule. The results of his studies indicate that 
transport of eroded material across the rings may be 
Significant for understanding the fine structure 
observed in the rings. 


ee, ore 983.  "Collisional Interaction of Ring 
Particles: the Ballistic Transport Process," Icarus 
eas 253-262. 


Ip develops a transport model independently of Durisen. 
The results of the two different approaches are quite 
compatible. Ip is able to assume partially elastic 
collisions, since he employs the Monte Carlo method in 
his simulation. He takes the collision probability to 
depend on optical thickness only; a more detailed 
See ree ote. Could casily™Se incorporated into his 
model, though. 


Seurces (Production Mechanisms) 


Mee lsom R.W. W980. "Photosputtering of Ice and Hydrogen 
Paowmma Saturn's Rings, Nature 283: 461. 


Soml=em argues that photodissociation of H90 On ring 
Bamueicles ~e@uld account for the relatively highspopu- 
lation of H atoms actually measured in the vicinity 
Simenews and B rings. Note that he calculates collision 
probability using l-e~*t with tT equal to one. He takes 
the sticking coefficient to be 0.22, but this value is 
based on laboratory experiments (Brackmann and Fite, 1961) 
MorcincO Mot Closely match the conditions in Saturn's 
rings. The overall influence of the uncertainty in 
these two factors should not affect the order of magni- 
~eucde Of Mas calculations, but probably limits the 
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precision of his estimation of total sroductlone. .e— | 
needed to produce the observed H atom population. 


Cheng, A.F. and Lanzerotti, L.J. 1978. £“IcCe Spe Gemeaae 
by Radiation Belt Protons and the Rings of Sa twine 
Uranus," Journal of Geophysical Research 83: 2597-2602. 


This paper is a detailed summary of sputtering effects 
ON Ging pareve wes: 


Cheng, A.F., Lanzerotti, L.J., and Pironelle,) Vee 
"Charged Particle Sputtering of Ice Surfaces in Saturn's 
Magnetosphere," Journal of Geophysical Research 87: 
4567-4570. 


Using information derived from the Pioneer II, Voyacemml 
and Voyager II missions, the authors reconsider the 
Sputtering mechanism, revising the result of their 
earlier (1978) paper. Since the actual ion fluxes 
observed were lower than those they assumed in 1978, 

the authors conclude that ion sputtering probably does 
not contribute significantly to the H atom atmoespiem 

tn the vicinity of the A and B rings of Sajeumne 


Ip, W-H. 1978. "On the Lyman-alpha Emission from the 
Vicinity of Saturn's Rings," Astronomy and Astrophysics 


TO eae 5 43 


This short paper suggests that H*™ ions escaping from 
Saturn's atmosphere may collide with ring particles 
where they are neutralized and re-emitted, thus 
contributing to the neutral H atmosphere in the vicinity 
Of “Ene  Fings-. 


Morfill, G.E. et al. 1983. “Some Consequences of 
Meteroid Impacts on Saturn's Rings," Icarus 55: 439=4a— 


Morfill et al. make careful estimates of the meteroid 

flux in the vicinity of Saturn. Using this esmimeeer 
coupled with a model for analyzing the effects of meteroid 
impact on ring particles, they proceed to make predictions 
about the ring halo, plamsa production, and the ring 
atmosphere. The authors conclude that meteroid impacts 
should contribute significantly to the ring) atmoscpaemes 


Ring Characetcerisredes 


i 


Cuzzi, J. et al. 1984. “Saturn's Rings: Properties vam 
Processes," in Planetary Rings, edited by Greenberg and 
Brahic, 75-199. University Of Arizenage «cee 


This paper is an up-to-date summary of information on 
Saturn's rings. It contd@ins a short Ssectronsde aed 
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to the ring atmosphere, which concisely explains 
current knowledge on that subject. Also, a special 
interest is an appendix which outlines radiative trans- 
fer theory (which is so important in deducing ring 
characteristics from observations). 


Pee est to, aw. eur ale 193842 >) "Saturn'’s Rings: Struc- 
ture, Dynamics and Particle Properties," in Saturn, 
edited by T. Gehrels and M. Matthews, 463-545. Univer- 
Sity of Arizona Press. 


This paper nicely complements the preceding paper 
KeuzZz4- et valee 1984) 


ee Marouf, E.A. et al. 1983. "Particle Size Distributions 
Maeoadturm Ss Rings from Voyager —T Radio Occultation, " 
Mearus 54: 189-211. 


This detailed paper provides the first precise estimates 

of particle size distributions ever given for the rings 

(results are calculated for four locations in the rings). 

The theory section in the paper also provides a clear 

application of radiative trasnfer theory. 

4. Weidenschilling, S. et al. 1984. "Ring Particles: 
Collisional Interaction and Physical Nature," in Planetary 
Rings, edited by Greenberg and Brahic, 367-415. Univer- 
Sity of Arizona Press. 


Besides providing an in-depth analysis of the effect 

of collisions between ring particles on the evolution 
of the rings, this paper gives information useful in 

the study of particle properties. 


Meeseeollaneous 


1. Chandrasekhar, S. 1960. Radiative Transfer. Dover 
Publications, Inc. 


Chandrasekhar produced the definitive work in the 

field of radiative transfer (radiative transfer theory 
seeks to specify the radiation field produced when 
incident radiation is scattered and absorbed by media 
consisting of discrete particles). This book is rigorous 
in approach, in fact it can be appreciated as a work in 
mathematical physics, alone. It is not light reading. 


Meeollack, J.B. 1975. "The Rings of Saturn,” Space 
Science Review 18: 3-93. 


Pollack's paper is a review of knowledge on Saturn's 
rings up to 1975. Much has been discovered since then, 
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rendering invalid some of the predictions in the paper 
This does not distract from the value of his work as 

a historical review. Reading this paper is an inspira- 
tional experience, because Pollack clearly demonstrag 
how knowledge about the rings increased through appli- 

cation of the scientific method. 


van de Hulst, H.C. 1981. Light Scattering iby een 


Particles. Dover Publications, Inc. 


This book is very readable. It essentially connects 
the mathematical physics of Chandrasekhar (1960) with 
the world of applied research. It provides the means 
by which the measurement of light scattering can be 
used to deduce the properties of the scattering medium. 
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